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The anharmonic response of charge-density wave (CDW) order to strong laser-pulse perturbations 
in lT-TaS2 and TbTea is investigated by means of a multiple-pump-pulse time-resolved femtosec- 
ond optical spectroscopy. We observe remarkable anharmonic effects hitherto undetected in the 
systems exhibiting collective charge ordering. The efficiency for additional excitation of the am- 
plitude mode by a laser pulse becomes periodically modulated after the mode is strongly excited 
into a coherently oscillating state. A similar effect is observed also for some other phonons, where 
the cross-modulation at the amplitude-mode frequency indicates anharmonic interaction of those 
phonons with the amplitude mode. By analyzing the observed phenomena in the framework of 
time-dependent Ginzburg-Landau theory we attribute the effects to the anharmonicity of the mode 
potentials inherent to the broken symmetry state of the CDW systems. 



Ultrashort laser pulses are a convenient tool for coher- 
ent excitation of phonons [l[ Q and collective electronic- 
lattice modes in CDW systemsH^. Due to availabil- 
ity of strong laser pulses the phonons can be driven 
far from equilibrium exposing anharmonic effects. The 
high excitation region has been already investigated 
with single-pump-pulse[3-[3 and double-pump-pulsefio'- 
12| sequences, mainly in elemental Bi, Sn and Te. 

In charge density wave (CDW) systems the phonon 
mode potentials are inherently anharmonic due to their 
coupling to the electron density modulation. [13] The an- 
harmonicity is the strongest for the Kohn-anomalv[l^ 
modes which become the amplitude and phase modes in 
the CDW state. While amplitude modes (AM) have been 
extensively investigated in the near equilibrium condi- 
tions by Raman[l4','l5^ and time resolved spectroscopy[3- 
0, [l^ and in the highly driven non-equilibrium condi- 
tions, where the CDW order is destroyeddli^, so far 
little attention[18] has been paid in the region inbetween. 

Here we report on an investigation of a new hitherto 
unexplored aspect of the CDW amplitude mode behavior 
under strongly driven non-equilibrium conditions in the 
ordered phase below the CDW photoinduced destruction 
threshold. Contrary to previous standard double-pump 
pulse (SDPP) high-excitation works floHl^ . where a pair 
of balancedfS^I pump pulses was used, we introduce a 
novel unbalanced double-pump-pulse (UDPP) approach 
in which we use the first and the strongest pump pulse 
(Pi) to excite large-amplitude coherent oscillations of the 
AM and other phonon modes and then use a standard 
pump-probe (P2-P3) pulse sequence to interrogate the 
system. By means of this novel approach we are able 
to directly investigate the anharmonicity of the effective 
AM potential, as well as detect the anharmonic coupling 
of the collective bosonic mode (the AM) of the CDW to 



other lattice modes. 

To establish generality, two layered chalcogenides 
which show different types of CDW ordering and also 
different electronic properties, were investigated: TbTca 
and ir-TaS2.[33il TbTca is a two-dimensional (2D) metal 
which shows an unidirectional incommensurate CDW 
state at the temperature used in our experiment (15 
K),[21, 22| while ir-TaS2 is in a commensurate insulat- 
ing CDW state at the relevant temperature (77K)[2i|. In 
both systems, in addition to the AM, several new Raman 
modes appear in the CDW state due to Brillouin-zone 
folding, [ill lillii Hi 



In our experiments the three pulse trains were derived 
from a 50-fs 250-kHz Ti:Al203 regenerative amplifier, 
with fvjj^ ~ 1.55 eV photon energy. The ps polarization 
was perpendicular to Pi and P2, which were parallel. To 
study the reflectivity change induced by the weaker P2 
pulse, AR2{t), we eliminate the Pi contribution, ARi{t), 
to the total transient reflectivity, AR{t), by means of the 
homodyne detection locked to the modulation of the P2 
pulse train. [The Pi pulse train was unmodulated as 
shown in FigfT^b).] 

In Fig. Wl^a.) we plot the raw UDPP photoinduced re- 
flectivity transients, AR2{t23)/R, in lT-TaS2 at differ- 
ent delays, ^12, between the pump pulses. The inten- 
sity of the P2 pulse train, I2, was set in the linear re- 
sponse region while the intensity of the Pi pulse train 
was 4-times larger corresponding to ^^30% of the CDW 
destruction treshold fluence. For comparison the raw to- 
tal photoinduced reflectivity transients AR{t23) / R, mea- 
sured in SDPP configuration, are shown in Fig. [IJc). 
In both cases the amplitude of the coherent oscillations 
periodically varies as ^12 is increased with a clear peri- 
odic suppression of the oscillations. Note that the sup- 
pression appears at different ti2 for each case. While 
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Figure 1: (Color online) AR/R transients as a function of ti2 
in ir-TaS2 in the UDPP configuration (a) as shown in the 
schematics (b), L - lens, S - sample, C - chopper, H.D. - ho- 
modyne detection system. For comparison AR/R transients 
in the SDPP configuration are shown in (c). UDPP power 
spectra in fr-TaS2 around the fundamental frequency of the 
strongest mode (d) and around the second harmonic of the 
fundamental frequency (e). The thin curve in (d) and (e) is 
the standard single pump-pulse spectrum. 



the linear SDPP effect [Fig. [^c)] is well knownfll, ^ 
[28| . and is understood as an interference due to the lin- 
ear superposition of two independently excited coherent 
oscillations [3, [23, [2^, the nonlinear effects observed by 
UDPP are completely new and hitherto undetected. 

In Fig. [TJd), (e) we plot the power spectra of the 
UDPP transients from Fig. [IJa). In addition to the AM 
mode with the frequency 2.41 THz we observe several 
weaker phonon modes above 3 THz and a weak second 
harmonic of the AM mode [see[IJe)]. The periodic inten- 
sity modulation of the modes, which strongly increases 
with increasing I2 (see Fig. [2]), is accompanied with a 
small periodic frequency shift. The shift is absent in 
the low-excitation SDPP configuration [see Fig. ^c) and 
(f)[. The modulation amplitude and phase vary among 
the modes and there is a ~ 7r/2 shift between the mod- 
ulation phases of the UDPP and SDPP case, similar to 
the observations in the high-excitation-density SDPP ex- 
periment in Te.(ll| In our case however, the phase shift 
persists down to the lowest excitation density. 



The periodic intensity modulation and the phase shift 
are observed also in TbTea, where in addition a beating 
in the ti2 dependence of the modulation amplitudes is 
observed (see Fig. [3^c)). 
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Figure 2: (Color online) Spectra of the strongest mode (a), 
(b) and weaker modes (d), (e) as functions of t\2 in the 
UDPP configuration at different intensities of the Pi pulse 
train in fr-TaS2. For comparison the low- excitation SDPP- 
configuration spectra are shown in (c) and (f). 
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Figure 3: (Color online) Integrated intensities of the strongest 
modes as functions of ti2 in UDPP configuration at Ji = 
4 X 72 in ir-TaS2 (a) and at h = 2.4 x h in TbTes (c). 
Normalized power spectra of the traces in the left panels (b) 
and (d). Open symbols correspond to the low- intensity SDPP 
response. Note the difference in the modulation frequency 
for the 3.38- THz mode in ir-TaS2 (b) and 2. 63- THz mode 
in TbTes (d) between the UDPP (full triangles) and SDPP 
(open triangles) cases. 

The modulation frequency of each mode [see Fig. [3] 
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(b) and (d)] is correlated to the respective mode eigen- 
frequency. Surprisingly, in the UDPP configuration there 
is also a clear cross-modulation of the 3.38-THz and 3.85- 
THz mode intensities with the 2.41-THz AM frequency 
and the 2.63-THz mode intensity with the 2.20-THz AM 
frequency in ir-TaS2 and TbTca, respectively. 

To understand the observed phenomena we start with 
the simplest Ginzburg-Landau expansion of the free en- 
ergy: 



F^Fn 



T 



l]\A\^ + -\A\* + g{t)\A\', (1) 



in terms of the normalized complex order parameter, A. 
Tc is the critical temperature and g{t) represents the ex- 
ternal laser excitation. Due to the symmetry g{t) can 
only couple to \A\'^. To describe the dynamics we intro- 
duce the T = AM frequency, wo, and the dimensionless 
damping, 7 — Awq/i^o, and obtain using p^: 



1 A+\AVA 



-g{t)A. (2) 



Since g{t)A can not excite phase fluctuations [3J|, only 
the AM needs to be considered. The dielectric constant 
depends in the lowest order on |ylp:|29|] 



e = eo + ci|A|^ 



(3) 



so the reflectivity change in the UDPP configuration is 
written as: 



Ai?2(i) = Ai?(t) - Ai?i(t) (X - \Ai{t)\ 



(4) 



where Ai (t) is the solution of equation ^ with excitation 
from the Pi pulse only and A{t) the solution with the 
excitation from both pump pulses. 

There are two terms in ([2]) that are of interest: g{t)A 
and The term g{t)A leads to a periodic modula- 

tion of the coupling of A to the P2 pulse when the oscil- 
lation amplitude after the Pi pulse is large. This effect 
however does not explain the main features of our obser- 
vations: (i) if the oscillation amplitude is large, a signif- 
icant AM-overtone intensity is expected due to \A\'^ in 
which is not observed in the experiment, (ii) because 
the laser excitation initially drives A towards zerojs^ an 
intermediate amplitude of the oscillations at ti2 = fol- 
lowed by a minimum at uJoti2 — 7r/2, where A(t) is closest 
to 0, is expected. Instead, the first minimum is observed 
around 0^0^12 ~ 37r/2 and the amplitude is the largest at 
ti2 = in lT-TaS2 and at woii2 — 7r/2 in TbTca. The 
anharmonic term in ^ remains therefore the only 

possible origin of the observed behavior. 

The solutions of equation ([2]) with 7 = can be rep- 
resented as closed periodic orbits in the phase space (see 
Fig. [3] (a)). Due to the anharmonicity the frequency of 
the orbit decreases with the oscillation amplitude. The 



P2 pulse transfers the system from the initial orbit, set 
by the Pi pulse, to a final orbit with a different frequency 
which depends on ^12. This results in a beating of the 
Ai?2(i) oscillations due to an interference of the second 
and the first term in the right hand side of ^ corre- 
sponding to the initial and final orbits, respectively. The 
frequency of the final orbit periodically oscillates with 
increasing ti2 [see Fig. |4] (b)]. Within a single initial 
orbit period there exist two delays, woii2 ~ '/'o +7r/2 and 
~ 00 + 37r/2, at which the only effect of the P2 pulse is 
a phase shift within the initial orbit. In vicinity of these 
delays the beating is very slow resulting in a slow in- 
crease of the Ai?2(i) oscillations amplitude with t. This 
slow rise is suppressed when 7 is finite since the oscil- 
lations die out before a significant phase shift between 
the orbits builds up and the Ai?2(i) oscillations ampli- 
tude remains small. As a result a periodic modulation 
of the Ai?2(i) oscillations amplitude is observed in the 
simulations with finite 7 [see Fig. 0] (c)]. In addition 
to the more intensive maximum, reminiscent to the ex- 
perimental observations, an additional weak maximum is 
observed within a single period corresponding to ^12 with 
the the larger final-orbit frequency [see Fig. |4](b)|. The 
weak maximum is strongly sensitive to the characteristic 
timescale of the external laser perturbation, rg,[36l| and 
is completely suppressed by setting Tg — 2tt/ujq as shown 
in Fig. HI (d). 

To improve agreement of the model with the exper- 
iment we extended the Ginzburg-Landau expansion ([!]) 
with the gradient term, ^^|9A/(92:|^,[23| to allow for space 
variations of the order parameter perpendicular to the 
sample surface due to the finite penetration depth. A, 
of the laser pulses. The result of a simulation with the 
inhomogeneous order parameter [Fig. |4] (e)] better re- 
produces the main experimental features for the AM. 
Similarly to the homogeneous case the weak maximum, 
which is not observed in the experiment, is absent only for 
long enough uJoTg 27r. In our experiment looTi ~ 0.27r, 
where rj is the length of the laser pulses, indicating that 
the coherent oscillations are not excited by the impulsive 
stimulated- Raman-scattering mechanism^[, but rather 
by the displacive one[lj, which can lead to Tg > ti. 

Counter intuitively, the experimentally observed small 
periodic frequency modulation of the AM is not repro- 
duced in any simulation with a finite damping so the 
origin of the effect is beyond the present model. 

Next we turn to analysis of the weaker modes with 
displacements denoted Qi. Similarly to the AM we 
can exclude the direct driving terms, \Qi\'^g{t) and 
di{QiA*)g{t), as the sources of the modulation due to the 
relatively small displacements. However, in addition to 
the AM also other modes are coupled to the CDW charge 
modulation. As a result, mode displacements Qi become 
mixed with the AM displacement [37| A and the effective 
mixed-mode potentials become anharmonic. This leads 
to the cross-modulation of mode intensities with the AM 
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Figure 4: (Color online) Orbits of eq. ([2]) in the phase space 
in the absence of damping (7 = 0) (a). The long arrow rep- 
resents the initial excitation by the Pi pulse. The short ar- 
rows represent the additional excitation by the P2 pulse which 
transfers the system to different orbits depending on t\2- Sim- 
ulated power spectrum as a function of t\2 in the absence of 
damping (b), with damping (7 = 0.01) and short excitation 
pulses, uJoTg = 0.27r (c), with damping and long excitation 
pulses uooTg = 2-n (d), with damping, long excitation pulses 
and finite laser penetration depth, A/^=16, (e). 



frequency in addition to the self-modulation at the par- 
ticular mode eigenfrequency. The relative amounts of the 
modulation at both frequencies depend on the couplings 
to the AM and to the external laser perturbation. If after 
the Pi pulse the amplitude of a particular mode is small 
the main contribution to the modulation is due to the 
cross-modulation at the AM frequency, as is the case for 
the 3.38-THz mode in lT-TaS2 and the 2.63-THz mode 
in TbTca. If, on the other hand, the amplitude is large 
enough a significant self-modulation is expected as in the 
case of the 3.85-THz in ir-TaS2. 

In conclusion, by means of a novel unbalanced double- 
pump-pulse time resolved optical spectroscopy, we were 
able to detect the inherent broken-symmetry-state an- 
harmonicity of the amplitude-mode effective potential in 
two distinct CDW systems. We also found a clear evi- 
dence of the anharmonic mixing of certain phonon modes 
with the amplitude mode originating from their mutual 
coupling to the electronic-density CDW modulation. We 
showed that the observed effects can be described in the 
framework of time dependent Ginzburg-Landau theory. 

This work has been supported by Slovenian Research 
Agency and the CENN Nanocenter. 
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